ttQZ JfcLL VOVN 




NATIONAL ADVISORY COMMITTEE 
FOR AERONAUTICS 

TECHNICAL NOTE 2044 

PRESSURE DISTRIBUTION AND SOME EFFECTS 

OF VISCOSITY ONJ3LE-NDER INCLINED 

BODIES OF REVOLUTION 

By H. Julian Allen 

Ames Aeronautical Laboratory 
Moffett Field, Calif. 



Washington 
March 1950 



NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


TECHNICAL NOTE 2044 


PRESSURE DISTRIBUTION AND SOME EFFECTS 
OF VISCOSITY ON SLENDER INCLINED 
BODIES OF REVOLUTION 
By H. Julian Allen 


SUMMARY 

In connection with a study of the flow around slender inclined hodie 
of revolution, a simplified solution for the pressure distribution on 
such a body has been developed. The solution should be suitable for 
bodies of high fineness ratio even at low supersonic speeds, provided the 
angle of attack is small. Comparison with low— speed experimental results 
indicates that the observed flow separation phenomena can be explained 
in the terms of the calculated pressure distributions together with the 
theory of oblique viscous flows. 


INTRODUCTION 

The longitudinal distribution of cross force on inclined bodies of 
revolution, which was primarily of interest to airship designers in the 
past, was solved simply and effectively by Max Munk (reference l). Munk 
showed that the cross force f per unit length on an arbitrary body of 
revolution in a nonviscous fluid stream was given approx ima tely by 

f sin 2a ' (1) 

I 

where >V 

q Q stream dynamic pressure 

dS/dx change in body cross-sectional area with longitudinal distance 
along the body 

a angle of inclination 

Tsien (reference 2) investigated the cross force on slender bodies 
of revolution at moderate supersonic speeds — a problem of more interest 
at the present to missile and supersonic aircraft designers — and showed 
that, to the order of the first power of the angle of inclination, the 
reduced Munk formula 
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f * 2<1 ° s a (2) 

was still applicable. This is not surprising when it is realized that 
the cross component of the flow field corresponds to a cross velocity 

V y Q “ V o sin a 

where V 0 is the stream velocity. Thus the cross component of velocity, 
and, hence, the cross Mach number will, for small angles of inclination, 
be a small subsonic value so that the cross flow will be essentially 
incompressible in character. 

Although the cross-force problem in a nonviscous flow has been satis- 
factorily investigated, the problem of determining the incremental pressure 
distribution due to inclined flow has not received such complete atten- • 
tion. Kaplan (reference 3 ) treated, in a thorough manner, the flow about 
slender inclined bodies, but the solution, which is expressed in Legendre 
polynomials, unfortunate ly , is tedious to evaluate. Laitone (reference 4), 
by linearizing the equations of motion, obtained a solution for the pres- 
sure coefficient on slender bodies of revolution in inclined flow, but, 
as will be seen, the solution, due to linearization, is inadequate in the 
general case. 

In this report, a simple incompressible potential flow method is 
developed for determining the incremental pressure distribution resulting 
from inclined flow on a slender body of revolution. The method, for the 
reasons discussed previously, should be applicable for subsonic and super- 
sonic flows wherein the cross velocity is small compared to the speed of 
sound. 


It. is known that viscosity actually plays an important part in the 
cross— flow phenomena even at moderate angles of attack. It may be shown 
theoretically that, under certain restrictions, the viscous cross force 
will be given by the drag of the body in cross-wise motion at a velocity 
equal to the crosswise component of the stream velocity. As a simple 
rule, the total cross force on a very slender body may be formed by adding 
this viscous cross force to the potential forces determined by Munk's 
formula. For such bodies this rule gives predicted characteristics in 
good agreement with experiment. A comparison is made in this report 
between the experimental pressure distributions on an inclined, model- 
airship hull and those calculated by the nonviscous theory. The compar- 
ison is instructive in indicating the conditions under which the calcu- 
lated characteristics will be obtained and in showing the manner in which 
viscosity effects are manifested. 
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NOTATION 

C constant of integration 

P surface pressure coefficient [ (p— P 0 )/q. 0 ] 

P a _Q surface pressure coefficient at zero angle of stream 
inclination [ (p a=0 -P o )A o ^ 

AP . incremental pressure coefficient due to angle of inclination 

[ (P-Pa^Ao 1 

f cross force per unit length on an inclined body of revolution 
L length of body 

M q free— stream ffach number 

p pressure on body or in field of body, as indicated 

p Q free— stream pressure 

P a _0 surface pressure at zero angle of stream inclination 
q. Q free— stream dynamic pressure 

r polar radius about body axis of revolution 

E radius at any station of body of revolution 
S cross-sectional area at any station of body of revolution 
t time 

Y q free— stream velocity 

V x axial velocity at any station of body surface 

V axial component of free-etream velocity 

o . 

cross component of free— stream velocity 
x axial coordinate 
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y ordinate in plane of inclination normal to axis of revolution 

z ordinate normal to plane of inclination and to axis of 

revolution 

a angle of inclination 

0 polar angle at out axis of revolution measured from the approach 

direction of the cross— stream velocity 

. -i cLR 
n tan — 

dx 

p fluid mass density 

cp velocity potential 


THEORY 

Consider the flow over the body of revolution shown in figure 1 which 


is inclined at an angle a 


to the stream of velocity 
V, 


V Q . If the body 


is slender, the axial component velocity V x at the body surface will not 

differ appreciably from the axial component V_ of the stream velocity. 

*0 

With this condition, it is clear that the flow may be treated by consider- 
ing the two-dimensional flow in a plane parallel to the yz plane, which 
plane is moving downstream with the constant velocity V x . In other 

words, the problem may be treated by determining the two-dimensional flow . 
about a circular cylinder which is first growing (over the forebody) and 
then collapsing (over the afterbody) vith t ime . 


The velocity potential for the cross flow at any 
in polar coordinates as 


x station is given 


q> = -V^ ( r + ^ ^ cos 0 (3) 

which in this moving reference plane is a function of time. 

Bernoulli’s equation for an incompressible flow which changes with 
time is 


p _ _ Sep _ 1 

p 3t 2 


(gj + (&J ] 


+ C 




Wow from equation ( 3 ) 
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tut 



<LR &R dz _ . 

— v tan Ti 

dt dx dt z o 


so that equation (5) becomes 


(5) 

( 6 ) 


It = ~^ Y Jo Y * 0 tan ( ?) 008 9 (7) 

Also, by differentiation of equation (3), 



so that equation (4) for the pressure at any point in the flow field 
becomes 


i = 2 Y J 0 \ tan ’I 


To 


■j^cos 2 0 1 


(I) 

- 071 


cos e — 


+ sin 2 9 


1 + 


(?/]■} 


+ C 


(9) 


for 


so 


r -> «, p -> p 


Po . V 


C = + 

P 2 

and hence equation (9) for the pressure at the surface of the body bee ones 
for r = R 


^2. = 2Vy o V Xo tan rj cos 9 + ^ 1 - 4 sin 2 0^ 


( 10 ) 
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and writing 


V = V sin a 
J o 


v x D = V o 008 a 


the surface pressure in coefficient form becomes 


p— p 

P = s 2 tan n cos 0 sin 2a + 

1o 


^1-4 sin 2 0 ^ 


4 sin 0 ) sin a ( 11 ) 


For small angles of inclination 

sin 2a ■ 2a 


sin 2 a m a 2 


then 


P = ^ 4 tan T) cos 0^ a + ^ 1 — 4 sin 2 0^ a 2 


( 12 ) 


This equation applies to a body of revolution of such extreme length 
that, at zero inclination, the local pressure at any x station is 
essentially the stream static pressure p Q . 


For bodies of lower fineness ratio at zero angle of inclination, 
the surface pressure at any station, -designated p a _Q, will differ 

slightly from the static pressure p 0 but, if the fineness ratio is not 
too low, the pressure p ass Q wil1 also closely approximate the pressure 
for some distance in this yz plane away from the body surface. Under 
this assumption that the pressure at the surface at zero inclination 
applies uniformly in the portion of the yz plane for which the major 
effects of the cross— flow distribution are felt, the change in pressure 
from p 0 to P a= o additive to, but will not otherwise influence, 

the cross— flow pressure distribution. Hence for any station on a body 
of high- f ine ness ratio for which, at zero inclination, the pressure is 
p a=0 , the pressure coefficient distribution at this same station under 

inclined flow conditions will be, from equation (ll). 


( 


P = Pq-o + V 2 tan ^ 008 0 


^ sin 2a + ^ 1 - 4 sin 2 0^ sin 2 a (13) 
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or, from equation (12) for small angles of inclination, 1 

P » P a _o + ^4 tan tj cos 0 ^ a + ^ 1 — 4 sin 2 0 ^ a 2 (14) 

The cross force per unit length of the body is then found as 

,2rt O Jt 


p n pit 

p cos 0 R (10 = 2q 0 / PR cos 0 d0 + 2 / 


p Q R cos 0 d0 


and clearly 


■r* 


R cos 0 d0 5 0 


Substituting P from equation (ll) gives 

n 


n n p* 

f = 2Rq Q P a _Q J cos 0 d0 + 4Rq Q tan tj sin 2a J cos 2 0 d0 + 

2Rq 0 sin 2 a 4 sin 2 0) cos 0 d0 

The first and third integrals are zero while the second integral yields 


and since 


then 


f = 2«Rq tan tj sin 2a 


2ttR tan tj = 2jtR — = — 
dx dx 


f = q ' — sin 2a 
° dx 


which is equation (l) derived by Munk for the cross force on slender 
airship hulls and, in the form. 




^Equation (l4), for the case in which tangent tj is constant, reduces 
to that derived by Busemann (reference 5) for the flow over an 
inclined cone. 
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that derived by Tsien for the cross force, to the order' of the first power 
of the angle of inclination, for slender bodies at moderate supersonic 
speeds. This development shows that these equations for the cross force 
are also correct to the second power of a. 


EXPERIMENTAL RESULTS AND DISCUSSION 
OF THE EFFECTS OF VISCOSITY - . 


In reference 6, a thorough Investigation at low speeds was made of 
the pressure distribution over a hull model of the rigid airship "Akron." 
Incremental pressure distributions due to inclination calculated by 
equation ( 13 ) for five stations along the hull at three angles of attack 
are compared with the experimental values in figures 2 to 6. In each of 
the figures Is shown a sketch of the airship which indicates the station 
at which the incremental pressure distributions apply. This comparison 
represents a severe test of the theoretical methods of this report since 
the method was developed on the assumption that the fineness ratio of the 
body is very large, while for the case considered the fineness ratio is 
only 5«9* 

At the more forward stations (figs. 2 to 4), the agreement is seen 
to be essentially good 2 but some discrepancy — particularly at values of 
9 near 180 — is evident which increases with increasing distance from 
the bow. The discrepancy increases rapidly, proceeding to sections 
(figs. 5 and 6) downstream of the maximum diameter section until (fig. 6) 
the entire distribution is affected. 

The disagreement that exists at the afterbody stations results from 
effects of viscosity not considered in the theory as will be seen from, 
the following: R. T. Jones, in reference 8, showed that for laminar flow 

on an infinitely long, inclined right circular cylinder, the behavior of 
the component flow of a viscous fluid in planes normal to the axis ^of 
revolution was independent of the component flow parallel to this axis. 
That is, viewed along the cylinder, the flow of a viscous fluid about the 
cylinder would appear identical to the flow about a section of a right 
circular cylinder in a stream moving at the velocity V Q sin a. Hence, 
separation of the flow would occur in the yz plane as a result of the 
adverse pressure gradients that exist across the cylinder. Jones demon- 
strated that this behavior well explained the cross forces on inclined 
cylinders that were experimentally observed in reference 9* .That such 
separation effects also occur on the inclined hull model of the "Akron" 
is also evident, particularly in figures 5 and 6. 

At stations extremely close to the bow the method must be inaccurate 
as evident from the work of Upson and Klikoff (reference 7 ). 
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While the treatment of reference 8 explains qualitatively the observed 
behavior of the flow field about the hull model considered, it cannot be 
used quantitatively,- in the general case, for at least three reasons: 

First, the theory of reference 8 was demonstrated only for laminar flows; 
the separation effects on bodies with turbulent flows may not' be the same 
as indicated in that reference. Second, the forward sections of bodies 
of finite length will behave more nearly as a circular cylinder set in 
motion initially from rest. Thus, although the adverse gradients exist, 
the flow, as shown in reference 10, will not have had sufficient time to 
exhibit the Usual separated flow characteristic of the steady— state flow 
across a circular cylinder. This, in part, explains the good agreement 
evident between the calculated and experimental values of figures 2 and 3. 
Third, the influence of the term 

2 tan q cos 0 sin 2a 

of equation (13) is to distort the typical circular cylinder pressure dis- 
tribution, given by the term . 

(1 — 4 sin 2 0) sin 2 a ' 

to move the calculated position of minimum pressure away from the 0=90° 
point and to change the magnitude of the pressure to be recovered. Over 
the forward stations of thie body where tan q is positive, the position 
of minimum pressure lies between 90° and l80° and the pressure required 
to be recovered is small and even zero, at the most forward stations. For 
the rearward station where tan q is negative, the minimum pressure lies 
between 0° and 90°, and the pressure recovery required is large and 
increases proceeding toward the stern. For the hull of the "Akron" model, 
the calculated line of minimum pressure along the hull is, for the angleB 
of attack of 6°, 12°, and 18°, as shown in figure J . 3 Since separation 
can only occur in an adverse gradient, it is clear that the line of 
separation will roughly follow the line of minimum pressures. Hence , 
again, the flow about forward stations will be, or will more nearly be, 
that calculated for a nonviscous fluid. Over the rearward stations the 
flow separation should tend to be even more pronounced than would occur 
on a right circular cylinder. 

From the foregoing, it is evident that the potential flow solution 
for the pressures on inclined bodies can only be expected to hold over 
the forebody,, and that over the afterbody the pressure distribution will 
be importantly influenced by the fluid viscosity. 

0 ' r r ( 

It is of interest to note in this figure that even for small angles of 
inclination the lines of minimum pressure become oriented close to 
the direction of the axis of revolution, while at zero inclination 
it must, of course, be normal to this axis. 
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The equations of this report for the pressure distribution on an 
inclined body were developed on the assumption that the flow was incom- 
pressible. Nevertheless, as noted earlier, the equations should be appli- 
cable to supersonic flow when the appropriate values for P a _Q are used 
as long as the cross component of the Mach number Mq sin a is not too 
large compared with the critical Mach number of a circular cylinder. 

In reference 4, the pressure distribution was determined for bodies 
of high fineness ratio to the order of first power in a and the results 
agree with equation (l4) of this report to this order. In light of the 
present solution of equation (l4), it" is clear that the solution of ref- 
erence 4 is only correct when a is not only small but small by compar- 
ison with q, for, when a is of the same order as q, the term 

(1 — 4 sin 2 0) a 2 

assumes equal importance with the linear term. Hence, for the usual small 
cone angles used on supersonic bodies, the solution of reference 4 is 
seriously limited in the range of angles of inclination for which it will 
give sufficiently accurate results. 


Ames Aeronautical laboratory. 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Dec. 28, 1949. 
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Figure 2. - Calculated and" experimental pressure distribution on 
a model hull of the U.S.S. Akron at station O.lll. 
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Station for incremental pressures, j- = 0.355 



— Theory 
o Experiment 


(a) oc. = 6 



(c) a = !8 . 
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8, deg 


Figure 4- Calculated and experimental pressure distribution on 
a model hull of the US.S. Akron at station -r= 0.355. 
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Figure 5- Calculated and experimental pressure distribution on 
a mode! hull of the U.S.S. Akron at station £ = 0.708. 
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Station for incremental pressures, = 0.867- 
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— Theory 

J^— o Experiment 


Separated flow 


(a) a = 6* 



(c) oc = 18° 
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Figure 6 -Calculated and experimental pressure distribution on 
a mode! hull of the U.S.S. Akron at station £=0.867 





